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ABSTRACT 

Recent progress in non-perturbative investigations of the electroweak phase 
transition is reviewed, with special emphasis on numerical simulations in the 
four-dimensional SU(2) Higgs model. 


1. Introduction 


The symmetry restoration in the electroweak Standard Model at high temper¬ 
ature! is an interesting feature in the history of the early Universe. The anomalous 
baryon number violating processes are fast for temperatures above the electroweak 
phase transition of symmetry restoration, and are getting very slow soon after it. It 
is also possible that the baryonic excess was produced during the phase transition 
itself!. The conditions for this are: a strong enough first order phase transition and 
a sufficiently strong CP-violation in the Standard Model. Concerning the strength 
of CP-violation, in spite of some optimism!, serious doubts have been raised!. In 
order to determine the properties of the electroweak phase transition one needs non- 
perturbative studies, which can be performed within the framework of lattice regular¬ 
ization. In this lecture a short review will be given about some recent developments 
in this field. The discussion will be centred around a recent series of numerical simu¬ 
lations performed in the four-dimensional SU(2) Higgs model! I. 


1.1. Dimensional reduction 

The main difficulty in numerical simulations of the thermodynamics of the 
electroweak plasma on the lattice has been identified in earlier studied!: as a conse¬ 
quence of the weak gauge coupling, the observed first order phase transition gets very 
weak with increasing Higgs boson mass ( Mh) already near Mh — Mw (with Mw the 
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IT-boson mass). The screening masses (inverse correlation lengths) near the phase 
transition become much smaller than the temperature. It becomes very difficult to 
deal with this two-scale problem in a numerical simulation. 

This suggests a reduction of the problem to a three-dimensional field theory 
by integrating out the heavy modes with non-zero Matsubara frequencies0~0. The 
simplest possibility is to consider the three-dimensional fundamental Higgs model 
with couplings determined by a perturbative matching procedure^. This amounts to 
choosing some perturbatively well calculable set of parameters of the four-dimensional 
continuum theory (usually derived from the effective potential in the phase with 
broken symmetry) and setting them equal to the corresponding parameters in the 
reduced three-dimensional continuum theory. Since the three-dimensional model is 
super-renormalizable, the correspondence of the renormalized parameters to the bare 
parameters in lattice regularization can be determined, near the continuum limit, in 
a two-loop perturbative calculationEl. In this way the non-perturbative properties of 
the phase transition can be numerically studied in lattice simulationsElo, which are 
less demanding than the direct lattice simulations in the original four-dimensional 
theory. 


Originally, the three-dimensional reduced model contains also a Higgs field in 
the adjoint (triplet) representation, but its mass is governed by chromo-electric screen¬ 
ing, and hence it is heavier than the Higgs doublet field and the three-dimensional 
gauge field, which gets its mass from chromo-magnetic screening. Therefore, at least 
in a first approximation, one can also neglect the adjoint Higgs field and stay with a 
simple three-dimensional fundamental SU(2) Higgs model. Further reduction, leading 
to scalar degrees of freedom alone, seems to give qualitatively different resultdil, and 
hence is not advisable. 

The obvious question about reduction is the effect of the necessary approxi¬ 
mations on the non-perturbative results. In order to discuss this, it is advantageous 
to think about the four-dimensional continuum theory as the continuum limit of a 
lattice-regularized model. Although the continuum limit of the SU(2) Higgs model 
(and also of the electroweak sector of the Standard Model) is believed to be trivial, 
the renormalized couplings go to zero only logarithmically with the lattice spacing. 
Therefore the continuum theory with non-zero renormalized couplings can be real¬ 
ized to a very good approximation along the lines of constant physics defined by the 
constant values of two dimensionless parameters: the mass ratio Rhw = Mh/Mw 
and the renormalized gauge coupling (& defined, for instance, from the screened 
(Yukawa) potential of external charges!^. The points of the lines of constant physics 
can be parametrized by the logarithm of the inverse IH-boson mass in lattice units: 
r = log(aMw)^ 1 - Due to the triviality of the continuum limit, every line of constant 
physics with non-zero couplings ends on the boundary of the bare parameter space 
at a finite value of r = r max , but if the Higgs boson mass is not close to the triviality 
upper limit0 ( M H ~ 9M W ), the value of r max can be very large. In such cases, since 
lattice artefacts go to zero by powers of the lattice spacing, the continuum theory is 
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very well approximated. 

In a given point of a line of constant physics, with large r, the heavy and 
high-momentum modes can be integrated out by consideringthe effective action of 
averaged fieldslH'IH or by applying block-spin transformation®. In fact, for the pur¬ 
pose of obtaining a three-dimensional lattice model it is simplest to take the block-spin 
approach and to perform so many block-spin transformations that in the imaginary 
time (or inverse temperature) direction only a single layer of points is left over. In this 
way a three-dimensional lattice action is obtained for the light and low-momentum 
modes, which is exactly equivalent to the original four-dimensional lattice action. 
(Here we assume the ideal situation of a perfect action, when no approximations are 
made in the blocking procedure.) 

The three-dimensional lattice action contains, in general, many couplings. Be¬ 
sides the simplest nearest neighbour couplings, there are also less local terms, with 
non-locality at least as large as the inverse temperature represented by the periodicity 
in the fourth dimension of the original four-dimensional lattice. How well this multi- 
parameter action can be approximated by the simplest local action corresponding 
to the discretization of the simplest local three-dimensional continuum theory is an 
open question. The connection between the multi-parameter action and the simple 
local action giving the best approximation will, in general, also depend on the lattice 
spacing. This can be seen, for instance, if the integration of heavv degrees of freedom 
is done in two-loop perturbation theory with momentum cut-ofST In fact, the con¬ 
tinuum limit in case of the three-dimensional block-spin action is not the same as the 
continuum limit for the simplest local three-dimensional action. Starting from dif¬ 
ferent points of some four-dimensional line of constant physics with increasing scale 
parameter r, the lattice spacing for the three-dimensional block-spin action is not 
changing at all. (The couplings in the action are, of course, changing.) One can also 
say that one should not take the continuum limit of an effective held theory beyond 
the scale where “new physics” appears. 

The danger in taking only the simplest local three-dimensional action is that 
the approximation given by it may be of different quality in different parts of the 
parameter space: the action optimized in the perturbative region, deeply inside the 
Higgs phase, can deviate more from the block-spin action near the phase transition 
and in the symmetric phase. In fact, it is plausible that in the symmetric phase, 
which is similar to the high-temperature phase of a non-Abelian pure gauge theory, 
the role of non-locality is more important. 

This shows that the numerical simulations of the electroweak phase transition 
in its unreduced four-dimensional form are important, even if finally the dimensional 
reduction, in some form, will turn out to be a good approximation. 

1.2. Lattice action 

The numerical simulations in the four-dimensional finite temperature theory 
have been restricted up to now to the SU(2) Higgs model. This is a good and suffi- 
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ciently simple theoretical laboratory, where the qualitative features of the symmetry 
restoring phase transition can be understood. The well known continuum action of 
this model is given in terms of the SU(2) gauge held W p {x) (fi = 1, 2, 3,4; r = 1, 2, 3) 
and the complex scalar doublet held <f> a (x) (a = 1, 2). The coupling parameters are: 
the (squared) gauge coupling g 2 and the quartic scalar self-coupling A. In addition, 
there is a single mass parameter belonging to the scalar held, which can be represented 
in the Higgs phase by the vacuum expectation value v. 

The Euclidean lattice action of the SU(2) Higgs model can be written as 



{ 2 Tr (PxPx) + A 


2 Tr ( ( Pt ( Px) ~ 1 


2 4 

~kJ2 Tt Ox+A UxpCPx) 

fi=l 


( 1 ) 


Here U Xfl denotes the SU(2) gauge link variable, U p i is the product of four I/’s around 
a plaquette, and <p x is a complex 2 ® 2 matrix in isospin space describing the Higgs 
scalar held. The bare parameters in the action are f3 = A/g 2 for the gauge coupling, 
A for the scalar quartic coupling and k for the scalar hopping parameter related to 
the bare mass square /Iq by (j,q = (1 — 2X)n~ 1 — 8. In what follows we set the lattice 
spacing to 1 (a = 1), therefore all the masses and correlation lengths, etc., will always 
be given in lattice units, unless otherwise stated. 

The large scale numerical simulations! - !, which will be discussed in the present 
lecture, have been performed mainly on the APE-Quadrics massive parallel computers 
at DESY, and in a smaller part also on the CRAY’s of HLRZ-Jiilich. The published 
results refer mainly to two Higgs boson masses: M H ~ 20 GeV and ~ 50 GeV, 
always assuming that M w = 80 GeV. Recently, an intermediate Higgs boson mass 
Mh — 35 GeV has also been considered!!, mainly in order to be able to compare the 
results with past and future results of numerical simulations in the reduced three- 
dimensional models. This work is still going on and will be published soonM 

In the optimized updating procedure for the creation of the sample of field 
configurations, besides the algorithms described in detail previously!, also scalar field 
reflections!! have been applied. In fact, this latter algorithm turned out to be very 
efficient in reducing the autocorrelations among field configurations, hence in reducing 
the necessary computer resources. 


2. Latent heat 

One of the most important physical quantities characterizing first order phase 
transitions is the latent heat, which is the discontinuity of the energy density Ae 
across the phase transition. In the SU(2) Higgs model on the lattice its dimensionless 
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ratio with the fourth power of the transition temperature T c can be obtained! from 





d\ d(3 ' 

— ■ A Q x - — ■ 6AP p i 


( 2 ) 


Here L t is the lattice extension in the inverse temperature direction, the partial deriva¬ 
tives of the bare couplings are taken along lines of constant physics, and 

L<P,xp = xTr {Px+pUx^x) i 


Qx — 


^Tr(< p+Vx) ~ 1 


P p i = 1 


-Tr U, 


pi 


( 3 ) 


The signs of the discontinuities such as Ae etc. are usually defined as differences of 
values in the symmetric phase minus the Higgs phase. 

For the numerical evaluation of the latent heat the knowledge of the lines of 
constant physics is required. Since in the Standard Model, for not very high Higgs 
boson masses, we are interested in weak couplings, estimates of A (r) and (3{t) can 
be obtained from the perturbative renormalization group equations. Similarly, k(t) 
can be estimated using one-loop lattice effective potentials®. These estimates can 
then be ckecked in numerical simulations by determining Rhw and g R numerically. 
If necessary, one can also further tune the bare parameters in order to reach constant 
values of Rhw an d g\. In practice, in the investigated parameter range, the pertur¬ 
bative estimates are usually sufficient. Assuming that the temperature is equal to 
T c , the scale changes are determined by the lattice extensions in the fourth direction: 
T c = 1 /L t . Taking L t = 2,3,4,..., we get for the changes in the scale parameter 
At = log(3/2), log(4/2),.... 

The other ingredients in Eq. are the discontinuities of the global quantities 
defined in Eq. [3]. In finite spatial volumes the discontinuities are round-off by inverse 
volume ((L x L y L z )~ l ) effects. These finite volume effects can be eliminated, in prin¬ 
ciple, by a cumbersome extrapolation to infinite volume. Another possibility is to 
use the metastability of the two phases in the vicinity of the phase transition. This 
implies that, in large enough volumes, during the Monte Carlo updating the held 
configuration remains in one metastable phase for a very long time. In this way one 
is able to determine expectation values characterizing single phases, and hence e. g. 
the required discontinuities. One has to have in mind, however, that in this way an 
independent precise determination of the location of the phase transition is necessary, 
because the metastability extends over a finite range of parameters. The uncertainty 
of the phase transition parameters appears in the errors of the discontinuities!!. 

First resultJll on Ae/T c 4 as a function of the Higgs boson mass are shown in 
Fig. 1, together with a comparison to two-loop resummed perturbation theory!!. As 
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can be seen, the results of the numerical simulation (shown by triangles) agree qual¬ 
itatively well with perturbation theory. The latent heat characterizing the strength 
of the first order phase transition is decreasing considerably between Mu ~ 20 GeV 
and M h ~ 50 GeV. 


3. Interface tension 


Besides the latent heat, another important physical quantity characterizing 
the first order phase transitions is the interface tension (a): the free energy per unit 
area of the wall separating the two phases. Many characteristics of the nucleation 
dynamics are determined by these two quantities. The radius of the critical droplets 
is proportional to cx/Ae. The nucleation rate is approximately determined by 


Rn = 


cr 


_ (.v/Tc? 

(Ae) 2 T c (A e/T^Y 


(4) 


The typical droplet distance is roughly proportional to Rn. In fact, the quantity R n 
defined in Eq. § can be taken as a good characterization of the strength of the first 
order phase transition: for R n = 0(1) we have a strong transition, whereas R n <C 1 
corresponds to a weak one. 

In numerical simulations one can exploit the tunnelling between the two nearly 
equal minima of the effective potential, which correspond to the two phases. On 
large enough lattices one can study the properties of the mixed phase when, as a 
result of tunnelling, both phases are simultaneously present in the system. For this 
purpose an elongated geometry of the lattice is advantageous. In a four-dimensional 
simulation, when the fourth dimension gives the inverse temperature (T = L 7 1 ), one 
can take an extension L z much larger than the other two “transverse” directions: 
L z L x = L y . In this case, for large enough transverse extensions when the mixed 
phase appears, there are transverse walls separating the two phases. For periodic 
boundary conditions, and not too large L z , there are just two walls and two regions 
(one for each phase). 

This situation can be used to determine the interface tension in several ways. 
For instance, one can determine a from the small energy difference in the transfer 
matrix between the asymmetric and symmetric combinations of the two phased®. 
Another possibility is to measure the doubly peaked distribution of some order pa¬ 
rameter, when a is given by the ratio of the maxima to the flat region between the 
two peaks. This can be effectively done by the “multicanonical” simulation method®. 
In the “two-coupling” method® one divides the long lattice extension into two halves 
with different bare parameters: in one half slightly above (k > k c ) in the other half 
slightly below (k < k c ) of the phase transition. In the range 20 GeV < Mh < 50 GeV 
all three methods are applicable and give compatible resultsi’i. In Fig. 2 the obtained 
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interface tension at Mh — 20 GeV and ~ 50 GeV is shown and compared with two- 
loop resummed perturbation theory®. At M fr ~ 35 GeV a recent result idl 


= 0.06 ±0.01. (5) 

35 GeV 

This shows that also the interface tension is a steeply decreasing function of the Higgs 
boson mass. 

The combination introduced in Eq. (4], which is a good characterization for the 
strength of the phase transition, is decreasing from R n ~ 0.23 at Mh ~ 20 GeV to 
R n ~ 0.00003 at M H ~ 50 GeV. That is, the phase transition at M H ~ 50 GeV 
becomes very weak indeed. 

The interface tension, more precisely cr/T Cl plays a decisive role also in deter¬ 
mining the minimal spatial lattice size required for an accurate numerical investiga¬ 
tion of first order phase transitions. In order to determine discontinuities of global 
averages, screening lengths in separated metastable phases and other similar char¬ 
acteristics of the transition, one has to reach a mixed phase situation, which can 
be characterized by a well developed double peak structure of the order parame¬ 
ter distributions. The double peak structure appears when L x L y a/T c = 0( 1), with 
L x , L y < L z the transverse lattice extensions. The lattice in this volume should be 
fine enough to also resolve the interface structure, which becomes thiner for increasing 
Mh . These are physical requirements setting a minimum on the spatial lattice size, 
which are practically more important than the conditions imposed by the smallness 
of O(a) lattice artefacts. This implies that the requirements on spatial lattice sizes 
are essentially the same in dimensionally reduced three-dimensional simulations as 
in unreduced four-dimensional ones. The gain in computational resources by using 
the simplest local three-dimensional action is essentially equal to a factor L t . The 
consequence is that, for instance, the extension of the present numerical studies to 
Higgs boson masses beyond Mh — 50 GeV is not much easier in three dimensions 
than in four. 


cr 

T> 

C 


4 . Discussion 

The large scale numerical simulations discussed in the previous sections show 
that the symmetry restoring electroweak phase transition is of first order, with fast 
decreasing strength, in the Higgs boson mass range 20 GeV < Mh < 50 GeV. 
This conclusion is in qualitative agreement with two-loop resummed perturbation 
theory®. At Mh — 50 GeV the strength of the phase transition is definitely not 
enough to produce a strong enough out-of-equilibrium situation in the early Universe 
for the creation of the observed baryon asymmetry. The observed difficulty of iden¬ 
tifying the Erst order phase transition near Mh — 80 GeV showJ^I that in the range 
50 GeV < Mh < 100 GeV presumably no surprises are to be expected. Nevertheless, 



before a definite conclusion, the non-perturbative studies have to be extended towards 
higher Higgs boson masses, which are actually the experimentally allowed ones in the 
minimal Standard Model. This seems feasible in both unreduced and dimensionally 
reduced SU(2) Higgs models. 
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